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ABSTRACT 

Cosmic-ray and X-ray ionisations in molecular gas produce a weak far-ultraviolet flux through the radiative 
decay of H2 molecules that have been excited by collisions with energetic electrons (the Prasad-Tarafdar mecha- 
nism). I consider the effect of this dissociating flux on the oxygen chemistry in C-type shocks. 

Typically a few percent of the water molecules produced within the shock front are dissociated before the gas 
has cooled to 50 K. The resulting column density of warm OH rises from 10^^ to 10^^ cm ^ as the ionisation 
rate is increased from 10^^^ s^^ (typical of dark clouds) to 10^^^ s^^ (adjacent to supernova remnants). These 
column densities produce substantial emission in the far-infrared rotational transitions of OH, and are consistent 
with the OH/H2O ratios inferred from ISO observations of emission from molecular shocks. For high ionisation 
rates the column of warm OH is sufficient to explain the OH(1720 MHz) masers that occur where molecular 
clouds are being shocked by supernova remnants. 

The predicted abundance of OH throughout the shock front will enable C-type shocks to be examined with high 
spectral resolution through radio observations of the four hyperfine ground state transitions of OH at 18cm and 
heterodyne measurements of emission in the FIR (e.g. from SOFIA) 

Subject headings: MHD — masers — molecular processes — shock waves — ISM: molecules — supernova 
remnants 
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1. INTRODUCTION 

The energetic electrons produced by cosmic-ray and X-ray 
ionisations in molecular clouds coUisionally excite the Ly- 
man and Werner bands of H2. The subsequent radiative de- 
excitations generate a weak flux of far-ultraviolet (FUV) pho- 
tons capable of dissociating many molecular species (Prasad & 
Tarafdar 1983). Models of chemistry in C-type shock waves 
have neglected the internally generated FUV photons, which 
are able to dissociate > 1% of each shock-produced molec- 
ular species before the gas cools to 50 K. Although this does 
not significantly modify the abundances of the parent species, 
the relatively small abundance of a dissociation product may 
be of interest. A particularly important example is the disso- 
ciation of water to form OH. It is well-established that H2O is 
formed efficiently from Ol in the hotter part of the shock front 
by endothermic reactions once T > 400 K (Draine, Roberge & 
Dalgarno 1983; Kaufman & Neufeld 1996). The dissociation 
of a percent of the water implies an OH abundance ^ 10^^, 
and column densities ^ 10^^ cm~^. 

This column of warm OH is sufficient to explain the 
high OH/H2O abundances inferred from ISO observations of 
shocked molecular gas associated with the HH54 outflow 
(Liseau et al. 1996), and the supernova remnant 3C391 (Reach 
& Rho 1998). Enhanced OH abundances associated with shock 
waves have also been detected at radio wavelengths. OH(1720 
MHz) masers are found where supernova remnants are running 
into molecular clouds (Frail, Goss & Slysh 1994; Yusef-Zadeh, 
Uchida & Roberts 1995; Frail et al. 1996; Green et al. 1997; 
Koralesky et al. 1998; Yusef-Zadeh et al. 1999a,b). Population 
inversion occurs by collisions in warm (40-125 K), moderately 
dense (^ 10^ cm~^) molecular gas in the absence of a signifi- 
cant FIR radiation field. In addition, a line-of-sight OH column 
between 10^^ and lO^^cm"^ is required for masing to occur 
(Elitzur 1976; Pavlakis & Kylafis 1986; Lockett, Gauthier & 
Elitzur 1999). The conditions for masing in this transition are 



most likely attained in the cooling tail of a C-type shock wave 
driven into the cloud by the overpressure within the remnant 
(Wardle, Yusef-Zadeh & Geballe 1998, 1999; Lockett et al. 
1999) provided that sufficient water can be dissociated in the 
shocked gas before it cools below ~ 50 K. 

In this Letter I present models of the oxygen chemistry in C- 
shocks that include dissociations by internally-generated FUV 
photons. In §H I outline the physical and chemical processes 
that are included in the calculations, and results for ionisation 
rates typical of dark clouds and clouds adjacent to supernova 
remnants are presented in §^ Warm OH column densities in the 
range 10^^-10^^ cm^^ are produced, consistent with the pres- 
ence of OH(1720 MHz) masers if the shock front harbouring 
the masers is propagating roughly perpendicular to the line of 
sight - as is expected because the masing column should have a 
small line-of-sight velocity gradient. The implications of these 
results are briefly discussed in §Q. 

2. SHOCK MODELS 

The shock structure is assumed to be steady and plane- 
parallel with the magnetic field perpendicular to the shock nor- 
mal. The medium is comprised of neutral and ionised fluids, 
with the magnetic field frozen into the ionised component and 
the ionised and neutral components coupled by elastic colli- 
sions. The ionisation fraction is low, so the inertia and thermal 
pressure of the ionised component are neglected. 

I follow Kaufman & Neufeld (1996) by writing the magni- 
tude of the drag force per unit volume on the neutrals due to the 
drift of the charged species as 
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where p and pi are the densities of the ion and neutral com- 
ponents, is the preshock value of pi, and Vd is the ion- 
neutral drift speed. The drag force is in the direction of the 
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ion drift through the neutrals, and there is an equal and opppo- 
site force per unit volume exerted on the ions by the neutrals. 
The coupling lengthscale L is determined by the detailed com- 
position of the ionised fluid, a mixture of molecular and metal 
ions, electrons and charged grains. The calculation of the up- 
stream ionisation balance by Kaufman & Neufeld (1996) yields 
L K. 10"-2(C/10""s"i)(nH/10'' cm"''') cm for the densities 
of interest here. 

The ionisation rate Q is generally assumed to be due to cos- 
mic rays, but X-rays effect molecular gas similarly (e.g. Mal- 
oney, Hollenbach & Tielens 1996) so both ionisation sources 
are notionally lumped together here. The cosmic -ray and X- 
ray ionisation rates depend on the environment: in dark clouds 
the standard interstellar cosmic-ray flux yields Q ^ 10"^^ s"^, 
whereas C ~ 10"^' s"^ in molecular gas adjacent to super- 
nova remnants. For example, for a remnant characterised by 
Lx = 10^^ ergs s^^ and a radius of 10 pc, C ~ 3 x lO-^^ s^^ 
(see fQ. The cosmic-ray ionisation rate may also be of this or- 
der if the low-energy cosmic ray flux is increased by a factor 
of 100 over the local interstellar value, as seems to be the case 
for GeV cosmic rays (Esposito et al. 1996). I therefore consider 
ionisation rates in the range 10~^^-10~^^ s"^. 

The ambipolar diffusion heating rate per unit volume, ne- 
glecting the heat capacity and radiative cooling by the ionised 
fluid, is Fvd', heating by PdV work, and ionisations are also in- 
cluded. Cooling is assumed to occur by collisional dissociation 
and vibrational transitions of H2 (Lepp & Shull 1982), and by 
rotational transitions of H2, H2O, and CO (Neufeld & Kaufman 
1993; Neufeld, Lepp & Melnick 1995). 

The shock models include a simple reaction network suffi- 
cient to follow the oxygen chemistry: 
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Rates for collisional dissociation of H2 were obtained from 
Lepp & Shull (1982); and vibrationally cold rate coefficients 
for (^-(||) were taken from Wagner & Graff (1987) apart from 
their forward rate for (^ which is ill-behaved below 300 K - 
the rate coefficient from the RATE95 database (Millar, Farqhar 
& Willacy 1997) is adopted instead. Photodissociation rates for 
H2O, OH and O2 by internally-generated FUV photons were 
obtained from Gredel et al. (1989). 

1 adopt shock parameters consistent with the OH(1720 MHz) 
observations towards W28, W44, and IC443, which imply 
post shock densities ^ 10^ cm"-^, and Zeeman field strengths 
~ 0.3 mG (Claussen et al 1997): a shock speed 25kms"^, 
preshock density nH= 10'' cm"-^ and preshock magnetic field 
strength 100 ^G. 

3. RESULTS 

The shock structure for an ionisation rate C, = 10"^^ s~^, 
typical of cosmic-ray ionisation in molecular clouds, is pre- 
sented in the top panel of Fig. |l|. The velocities are plot- 
ted in the shock frame, in which the shock is stationary, with 



the unshocked gas flowing in from z < at the shock speed 
(i.e. 25kms~^, being decelerated within the shock front and 
departing downstream towards z > 0. The middle panel 
shows the chemical abundances obtained when dissociations 
by internally-generated FUV photons are neglected. As has 
been found previously (Draine et al 1983; Kaufman & Neufeld 
1996) the preshock atomic oxygen is almost entirely incorpo- 
rated into water, with OH appearing briefly as an intermediate 
step. For comparison, the lower panel shows the effect of in- 
cluding dissociations by internally-generated FUV. The overall 
chemistry is not drastically affected, but dissociation of water 
produces an increasing abundance of OH downstream, with a 
significant column having accumulated by the time the gas has 
cooled to 50 K. Note also the gradual reappearance of Ol down- 
stream through the dissociation of OH and O2. 

The effect of raising the ionisation rate to 3 x 10"^® s"^, as 
produced by X-rays in a molecular cloud adjacent to a SNR, is 
illustrated in Fig. |[ The ionisation fraction of the preshock gas 
is increased by a factor of V30, reducing the thickness of the 
shock transition and the time scale for gas to flow through the 
shock front by the same factor The peak neutral temperature 
increases because the energy dissipated in the shock must be 
radiated away on a shorter time scale. Heating associated with 
ionisations, which is unimportant within the shock front, in- 
creases 30-fold and so the final postshock temperature increases 
slightly. The photodissociation rate is increased 30-fold, with a 
concomitant increase in the production of OH and Ol. 

Fig. 1^ summarises the production of OH within the shock 
models by plotting the run of OH column density with temper- 
ature within the shock front for different ionisation rates. An 
initial 'burst' of OH column is produced as oxygen is rapidly 
converted to water in the high-temperature portion of the shock 
front. The column produced during this phase decreases as the 
ionisation rate is increased because the OH destroying forward 
reaction in (Q) is strongly temperature dependent, and the tem- 
perature within the shock rises more rapidly for high ionisation 
rates because of the increase in the rate of ambipolar diffusion 
heating. If the FUV dissociation is neglected no more OH is 
produced (dashed curves); otherwise OH starts to be produced 
once the temperature of the shocked gas has dropped suffi- 
ciently so that the forward reaction in (Q) has slowed sufficiently 
to permit photodissociations to become effective. For high ion- 
isation rates the increase in the photodissociation rate (oc C) is 
partially offset by the decrease in flow time scale (oc C~^^^); 
thus NoH C^/^ in the cooling gas. 

4. DISCUSSION 

The line-of-sight column density of warm OH through a C- 
type shock is A^oh sec 6 where 9 is the angle between the shock 
normal and the line of sight. For ( > 10"^® s"^ and 9 ^ 80, 
NoH is sufficient to permit the formation of OH(1720 MHz) 
masers, thus typical X-ray fluxes near SNRs are able to produce 
the warm OH column necessary for OH(1720 MHz) masers. 

Lockett et al. (1999) have argued that the X-ray flux from 
SNRs is < 0.02 of the level that could produce the OH col- 
umn density necessary for masing in the 1720 MHz line. This 
discrepancy can be traced to the conversion from X-ray en- 
ergy flux Fx to ionisation rate. Both Wardle et al. (1998) 
and Lockett et al. (1999) used the conversion C(s~^) « 3 x 
10~^''Fx(ergcm"^ s"^) (Maloney et al. 1996), but this as- 
sumes a hard X-ray spectrum appropriate to AGNs. The softer 
(E ~ IkeV) spectrum from the hot gas in SNRs produces 
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Fig. 
nil = 



z (10^® cm) 

1. — (Top panel: velocity and temperature profile (the temperature scale is on the right) of a 25kms~^ C-type shock propagating into gas with 



10^ cm ^ and iJo = lOO/i G. The X-ray + cosmic ray ionisation rate is assumed to be 



10- 



s ^ (see text). Middle panel: oxygen chemistry 



in the absence of dissociation by the FUV radiation field induced by ionisations. Lower panel: oxygen chemistry with dissociations included. 
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Fig. 3. — Solid: OH column density versus temperature within C-type shock models with = 25kms~^,nH = 10'' cm"'', and Bo = 100/^G. The curves 
are labelled by log s~'). Dashed: The same, but with photodissociations neglected. The dotted box indicates the conditions necessary for 1720 MHz masers to 
form assuming that the shock is propagating at 80° to the line of sight. 
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many more ionisations per unit X-ray luminosity, as ionisa- 
tions are dominated by low-energy X-ray photons. In this 
case, ( ~ NecrFx where Ne ~ 30keV~^ is the mean num- 
ber of primary and secondary electrons produced by the ab- 
sorption of unit energy, cr w 2.6 x 10^^^ cm^^ is the pho- 
toabsorption cross-section per hydrogen nucleus at 1 keV, and 
Fx = Ljf/47ri?2. With Lx = lO^^ergs^i and R = 10 pc, 
this yields <^ ^ 4.6 x 10~^^ s^^ provided that the hydrogen col- 
umn density Nh < 10^^ cm^^. Note that externally-incident 
UV may make a similar contribution to the OH column if Ay 
is not too large (Lockett et al. 1999). 

The enhancement of OH described in this Letter permits the 
structure and kinematics of C-type shock waves to be stud- 
ied through the modelling of the excitation of OH throughout 
the shock front, and the calculation of emission and absorp- 
tion line profiles. Supernova remnants with associated 1720 
MHz masers (Frail et al. 1994, 1996; Yusef-Zadeh et al. 1995, 
1996, 1999a,b; Green et al. 1997; Koralesky et al. 1998) are 



obvious observational targets. At radio wavelengths, the four 
ground-state transitions at 1612, 1665, 1667 and 1720 MHz 
can be observed with sub-km/s velocity resolution in absorp- 
tion against the background continuum from the remnant. For 
example, W28 shows extended OH absorption around the OH 
maser positions (Pastchenko & Slysh 1974; Claussen et al. 
1997). The predicted warm OH column within C-type shocks 
implies substantial far-infrared emission in low-lying OH rota- 
tional transitions. This will, for example, be easily detectable 
by the GREAT heterodyne spectrometer planned for SOFIA, 
an instrument that will be capable of 0. 1 km/s or better velocity 
resolution. 

I thank Cecilia Ceccarelli, Ewine van Dishoeck, Michael 
Kaufman and Mark Wolfire for discussions, and Phil Lockett 
for comments on the manuscript. The Special Research Centre 
for Theoretical Astrophysics is funded by the Australian Re- 
search Council under its Special Research Centres programme. 
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